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Abstract 

We compute the critical temperature T c for the deconfinement 
transition of pure QCD on coarse lattices, with Nt = 2, 3, 4, and lattice 
spacings from .33 fm to .15 fm. We employ a perturbatively improved 
gluon action designed to remove order a 2 and a s a 2 errors. We find 
that T c in units of the charmonium IP-IS splitting and the torelon 
mass is independent of a to within approximately 5%. 



PACS number(s): 12.38.Gc, 11.10. Wx, 12.38.Mh, ll.10.Gh, 05.70.Fh 



1 Introduction 



Recent calculations using improved lattice actions have demonstrated in sev- 
eral systems that accuracies of several per cent could be achieved with lattice 
spacings a as large as half a fermi. Use of an improved gluon action designed 
to remove order a 2 and a s a 2 errors demonstrated rotational invariance for 
the static quark potential independent of spacing errors to within 5% on lat- 
tices with a = .4 fmQ. Heavy quarkonium simulations with a nonrelativis- 
tic action corrected for order v 2 relativistic effects and a 2 errors produced 
masses with systematic errors as low as 5 MeV [@, [l[]. Other applications 
have included glueballs||, B mesons [Q, light fermions 0, || and QCD at 
high temperatures 0. 

The cost of lattice simulations typically increase as 1/a 6 ; improving accu- 
racy by decreasing a rapidly becomes prohibitively expensive. An attractive 
alternative is to add effective interactions to remove spacing errors directly 
in the lattice action, bringing it closer to the continuum || . The use of 
tadpole- improvement for the link variables and renormalized lattice pertur- 
bation theory in the design of these interactions are essential ingredients in 
these recent successes 0. Because these corrections can be computed per- 
turbatively, parameters need not be tuned, nor new ones introduced beyond 
those of continuum QCD |L0|. 



The efficiency with which improved actions remove lattice errors dra- 
matically increases the range of calculations accessible to currently available 
resources. In this work we employ an improved gluon action to compute 
the deconfmement temperature of QCD. In addition to its intrinsic impor- 
tance as a fundamental nonperturbative prediction, it provides an excellent 
quantity with which to test the accuracy of this and other improvement 



schemes 11, 12, 13 



2 Improved Gluon Action 

We used the improved action of Ref . |L4], |l|, []J , 

3[U] = /3 p i£§ReTr(l-[/ pl ) 
pi 

+ /3 rt £§ReTr(l-[/ rt )+/3 pg £§ReTr(l-tf pg ), (1) 

rt pg 



2 



with U p \ the Wilson loop about a single plaquette, U It a loop around a 1 x 2 
rectangle, and U pg a parallelogram with links running along the opposing 
edges of a cube. These terms are sufficient to provide the dimension four and 
six operators needed for an effective action accurate to order a 4 . 
The coefficient (3 p \ is the only input parameter, while 

&t = -7^2 (! + °- 4805 ««) > A>g = °- 03325 ( 2 ) 

were computed in terms of /3 p i using tadpole-improved perturbation the- 
ory |16| . The division of each link by uq gives tadpole-improved operators. 
These ensure that the connection between these quantities and their con- 
tinuum equivalents are not ruined by large tadpole contributions, and are 
essential for the effectiveness of lattice perturbation theory and the improved 
action |§. The measured expectation value of the plaquette determines both 
Mo and the renormalized coupling constant a s self-consistently, with 

U „=(iReT r(C / pl )) 1/4 , a , = _Ki5f^M. (3) 
V3 \ wj ' 3.06839 v ; 



3 Deconfinement Transition and 
Critical Temperature 

Green's functions computed on a Euclidean lattice periodic in time with finite 
extent t are equivalent to those at the finite temperature T — 1/t, averaged 
with the usual Boltzmann weight exp(—H/T). In particular, the correlation 
of a Polyakov loop or Wilson line, 

P(x)=TrntW~, (4) 

t 

which wraps around the lattice in the time direction, and its adjoint are 
related to the free energy Fqq of a static quark- ant iquark pair separated by 
x according to 

(P(x)pt(o)) oc exp(-F Q Q/T). (5) 
At large separation, cluster decomposition requires 

(P(x)pt( ))^|(P(0))| 2 . (6) 
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In the confining phase, Fqq increases with |x| and (-P(O)) must therefore 
vanish; in the deconfined phase, P may acquire a nonvanishing expectation 
value. Pure SU(3) lattice gauge theory, including improved actions, possesses 
a global symmetry under multiplication of timelike links at a single time slice 
by the elements of the group center, Z 3 . Nonzero (P) indicates the breakdown 
of this symmetry. Therefore (P) is a suitable and relatively simple order 
parameter with which to study this transition 17]. 

As the temperature exceeds its critical value T c , the system moves from 
confined to deconfined phases. Measurements of P, which were clustering 
around zero, now cluster about the real axis, or to the two axes along the 
other Z3 elements, exp(±i2-7r/3). During the transition, the system will fluc- 
tuate between phases. We determine T c by the location of the peak in the 
susceptibility 

X = iV 3 «P 2 ) - <P> 2 ) (7) 

which measures the deviation in (P) these fluctuations induce. Here P is 
averaged over spatial lattice sites. 

Because the volume is finite, the transition temperature is smeared, lim- 
iting the precision with which we may determine T c . Also, above T c , given 
sufficient simulation time the system will tunnel between the three degener- 
ate vacua, and (P) will still average to zero. It is therefore useful in practice 
to fold these three axes together [ 18], |TT|1 , replacing P in Eq. (|7|) with Q, such 
that 

( ReP ; arg P e (— tt/3, tt/3) 

n=l PePe i2 ^ 3 ; arg P e (vr,-vr/3) (8) 
[ Re P e - i2 ^ 3 ■ argPe (tt/3, tt) . 

We move our lattice system through the phase transition by holding N t = 
1/aT fixed while adjusting /3 pl , and therefore implicitly the lattice spacing a. 
We are measuring a quantity defined at infinite spatial volume, and N t should 
be small relative to the number of spatial sites per side, N s . We determined 
(3 c pl on lattices with (iV 3 x N t ) of (4 3 x 2), (6 3 x 2), (6 3 x 3), (9 3 x 3) and 
(8 3 x 4). 

After determining the critical coupling (3^ at which the transition occurs, 
we computed a convenient mass M with this same but at zero temper- 
ature, with N t > N s . The value of Ma from this simulation fixes a, and 
therefore T c , in units of M. Comparison with data for M gives T c in physi- 
cal units. We chose for M the charmonium IP — 15 splitting, and also the 
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torelon mass, which is related to the string tension. 

To determine we began with (3 V \ below transition and used one thou- 
sand Metropolis updates to allow the lattice to thermalize. We then sampled 
the Polyakov loop after every five updates for the next one or two thousand. 
After stepping (5 V \ up and allowing several hundred more updates, we began 
sampling again, repeating until (3 V \ had crossed its critical value. To ensure 
that the number of updates was sufficient to thermalize, we stepped f3 v \ back 
down through the transition to check for hysteresis. 

Figs. HI through |] display the susceptibility plots for each lattice as a 
function of f3 v \. The up and down arrows indicate data from upward and 
downward sweeps. The solid and dashed lines show independent fits of each 
to a Gaussian and a first order polynomial. The error bars indicate bootstrap 
estimates. We determine (3^ from the location of the peak, and assign the 
Gaussian width as the error. Table |T] presents the weighted average of (3^ for 
each T c . 




6.65 6.75 6.85 



7.15 7.25 




6.85 6.9C 



i.95 7.00 7.05 



Figure 1: Susceptibility x f° r 4 3 x 2 Figure 2: x f° r 6 3 x 2 lattice, 

lattice. Arrows indicate upward 
or downward sweeps of f3 p \ through 
transition. 

While we found locating the peak in x the most reliable method for 
locating an alternative is to locate the phase change by directly observing 
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Figure 3: x f° r 6 3 x 3 lattice. 



Figure 4: x f° r 9 3 x 3 lattice. 
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Figure 5: x f° r 8 3 x 4 lattice. 



pi 



1/2 1/3 1/4 



6.99(1) 7.49(1) 7.78(2) 



Table 1: Critical coupling (3^ as a function of critical temperature. T c is in 
lattice units. 



P as a function of (3 P \. Ref. [I9| implements this by choosing a threshold 
in the fraction of (-P)'s found near any Z3 axis. Fig. |6] displays the real 
and imaginary parts of P measured during an upward sweep of /3 p i on the 
9 3 x 3 lattice. The transition is evident. Below (3 P \ of about 7.45, P clusters 
about the origin. During the transition it explores the three degenerate vacua 
which break Z 3 , before settling into one above the transition. Given sufficient 
simulation time, P will populate each arm equally. In our study, it rarely 
explored more than two of the three. 



4 Charmonium 

To test the dependence of these results on the lattice spacing and to convert 
the critical temperatures to physical units, we computed the mass difference 
Mps between the IP and IS states of charmonium on lattices with the same 
spacings at zero temperature. We chose this quantity because it is experi- 
mentally well-measured, simple to compute accurately using a nonrelativistic 
action (NRQCD) |20], [§, small and therefore insensitive to finite lattice er- 
rors, and its splittings are essentially independent of the charm quark mass. 
We used an NRQCD quark action which removes the leading order a 2 and 
v 2 errors in lattice spacing and quark velocity, and so is consistent with the 
gluon action. The dominant uncorrected error is 0(a 2 v 2 ), due to the a 2 error 
in the lattice version of the D • E relativistic 0(v 2 ) correction. We did not 
include dynamical light quarks. Spin splittings are absent to this order, and 
we compare our results to spin-averaged data. 

Table |2] presents the values for Mpg for these lattices. We tuned the charm 
quark mass Mq by requiring that the energy for states carrying nonzero 
momentum satisfy the usual nonrelativistic dispersion relation. Comparing 
the lattice result to the experimental spin-averaged value for this splitting, 
Mpg = 457.9 MeV, fixes the inverse lattice spacing. 

We obtained these masses by 2 x 2 correlated matrix fits of exponentials 
to three propagators for which either the source, the source and sink, or 
neither were smeared by the lattice equivalent of exp(— p 2 /12) ~ (1 + V 2 /12), 
designed to suppress states with large momentum. We fit to both the ground 
and first excited states in each channel. Including the excited state prevented 
it and higher states from contaminating the value for the ground state, and 
allowed us to include data from lower times. We found this procedure useful 
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T c 


lattice 


M c 


M PS 


a- 1 [GeV] 


6.99 


1/2 


3 3 x 12 


2.0 


.77(1) 


.596(8) 


6.99 


1/2 


4 3 x 12 


2.0 


.794(4) 


.577(3) 


7.49 


1/3 


6 3 x 12 


1.3 


.494(6) 


.93(1) 


7.78 


1/4 


6 3 x 12 


.98 


.378(7) 


1.21(2) 


7.78 


1/4 


8 3 x 12 


.98 


.375(4) 


1.22(1) 



Table 2: Mps as a function of (3 pl . Errors are statistical. The critical temper- 
ature T c , the charm mass Mc, and the IP — IS splitting Mps are in lattice 
units. The experimental value for Mpg fixes a. 



P C P 1 


T c 


lattice 


T C /M PS 


T c [MeV] 


6.99 


1/2 


3 3 


x 12 


.651(14) 


298(6) 


6.99 


1/2 


4 3 


x 12 


.630(10) 


288(5) 


7.49 


1/3 


6 3 


x 12 


.675(13) 


309(6) 


7.78 


1/4 


6 3 


x 12 


.661(17) 


303(8) 


7.78 


1/4 


8 3 


x 12 


.667(14) 


305(7) 



Table 3: Critical temperature T c in units of Mpg and in MeV. 

in Ref. [fj], and especially helpful in extraction of torelon masses, discussed 
below. Typically we obtained good fits by including minimum times from 3 
to 6, to maximum times from 9 to 12. 

Table |3] displays the critical temperature in units of Mps and in MeV, 
using the experimental value for Mps- Fig. |7| presents the average T c for 
each a, around the central value T c = 300(6) MeV. Along with the statis- 
tical errors of quoted for M P g in Table §, we include errors associated with 
uncertainties in the critical value for (3 P \ and in the value for the average 
plaquette u , which is used to tadpole- improve the operators in the action 
which remove the 0(a 2 ) errors. To determine the former, we used one-loop 
perturbative evolution to estimate the variation of Mps with /3 P \. For the 
latter, a variation in uq would affect the ability of the improved action to 
remove 0(a 2 ) errors. We estimated the percentage error by 8uq/uq times the 
typical (and similar) 0(a 2 ) errors for these systems: a 2 (p 2 ) for charmonium, 
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and a 2 a 2 for the torelon, discussed below. Here a is the string tension, and 
5uo the difference between the input and measured value. We found these 
methods to provide a reasonably reliable estimate of the variation of Mpg 
and the torelon mass in the data over a fairly broad range in (3 V \ and uq. 

The values for T c are nearly consistent within these errors, and show 
little dependence on lattice spacing. We note that values for T C /M PS for the 
coarsest lattices (/^ = 6.99) deviate by about two standard deviations. This 
is a possible indication of the 0(a 4 ) systematic errors for which the action 
is uncorrected. These might be expected to range from about 5-10% for the 
coarsest lattice (about 1/3 fm) to a couple per cent for the finest (1/6 fm). 



3CC 

T c (MeV) 



Figure 7: Critical temperature T c as a function of lattice spacing a, around 
central value. Mpg sets the physical scales. 



5 Torelon Mass 

As a second test of scaling, we measured the torelon mass \x for torelons of 
two different physical lengths at these same (3 V \. The torelon is created by a 
Wilson loop strung along the full spatial extent of the periodic lattice, and 
averaged over spatial orientations and transverse positions [^l]]. As such, it 
is sensitive to both the string tension a and lattice length L. 
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/ft 


T c 


lattice 




ft* 


6.99 


1/2 


3 3 x 12 


1.46(2) 




6.99 


1/2 


4 3 x 12 


2.26(5) 




7.49 


1/3 


6 3 x 12 


1.42(1) 


3.0(4) 


7.78 


1/4 


6 3 x 12 


.676(6) 


1.8(3) 


7.78 


1/4 


8 3 x 12 


1.05(2) 


2.1(1) 



Table 4: Torelon mass jj and excited torelon mass jj* vs. /ft. Errors are 
statistical. T c is the critical temperature in lattice units. 



We obtained the torelon mass from the exponential decay in the correla- 
tion between two such loops at different times. More precisely, we accounted 
for periodicity in time by fitting to hyperbolic cosines rather than exponen- 
tials, and typically fit to time steps from one to six. For sufficient length L, 
the mass should go as 

H ~ oL , (9) 

reflecting the cost in energy of keeping static charges separated by this dis- 
tance. We will use /i to determine a value for a below. 

We may estimate the energy expected for excited states by comparing to 
those of a classical vibrating periodic string, 2iin/L. This would suggest first 
excited states split from the ground states by 2ti/N s in lattice units, close to 
the values we observe. For larger lattices, it therefore becomes increasingly 
important to include these states in fits. 

To extract the torelon masses, we performed correlated 3x3 matrix fits 
of three hyperbolic cosines to six torelon propagators which were smeared 
once, twice or unsmeared at either the source or sink. To smear these, we 
used the operator exp(— apj_) on the link variables ]22| , which suppressed 
contributions with large momentum transverse to the loop. Specifically, we 
applied n times the lattice version of (1 + aVj_/n), with n typically 7 and 
14, and a effectively .6 and 1.2. We reuniterized the link variables between 
applications. 

Table [| lists the results of these fits for the ground state masses and 
several of the excited states. For smaller lattices, we were unable to get 
accurate fits for the excited states, presumably due to the larger splitting, 
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6.99 


1/2 


3 


.343(7) 




1.79(4) 




7.78 


1/4 


6 


.370(7) 


.139(23) 


1.90(4) 


4.8(8) 


6.99 


1/2 


4 


.221(6) 




2.78(8) 




7.49 


1/3 


6 


.235(4) 


.111(14) 


2.87(4) 


6.1(8) 


7.78 


1/4 


8 


.238(6) 


.119(6) 


2.85(6) 


5.6(3) 



Table 5: Critical temperature T c in units of the torelon masses /i and /!*, 
and the ratio of these to Mps. Entries are grouped according to torelons of 
the same physical length L, given in lattice units. 



p c pl 


T c 


L 


L [fm] 


era 2 




T c /v^ 


^ [MeV] 


6.99 


1/2 


3 


.99(2) 


.603(1) 


2.33(2) 


.644(5) 


463(8) 


7.78 


1/4 


6 


.98(3) 


.142(2) 


2.26(2) 


.664(5) 


456(9) 


6.99 


1/2 


4 


1.37(2) 


.631(15) 


3.18(4) 


.630(7) 


458(8) 


7.49 


1/3 


6 


1.28(2) 


.266(4) 


3.09(2) 


.647(5) 


478(6) 


7.78 


1/4 


8 


1.29(3) 


.148(4) 


3.07(4) 


.651(8) 


469(6) 



Table 6: String tension a defined by Eq. (|i~0|), given in the dimensionless 
combinations with T c and L, and in MeV. Mps fixes the physical scales. 
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and settled for 2x2 matrix fits to two hyperbolic cosines. Our ground state 



masses fi fell below those of Ref. [12], from 3% below at T c = 1/2, increasing 
to 12% below for T c = 1/4. This may indicate excited state contamination 
in their fits, which would account for the increasing discrepancy with torelon 
length. We did reproduce their result for the torelon on a 3 3 x 16 lattice 
using an unimproved Wilson action at (3 — 5.10, obtaining 1.82(3), which we 
computed as a check. In this case we were unable to reliably fit the excited 
state or include time steps much beyond the first couple. 

In Table [5] we present the critical temperature in units of the torelon 
masses, as well as the ratio of these masses to Mps- The torelons we studied 
had two different physical lengths, determined by the ratio of L to 1/T C , and 
are grouped accordingly. In most cases, these quantities are independent 
of a to within statistical errors. There is roughly a 5% discrepancy for the 
coarsest lattice, which may indicate an 0(a A ) effect. Ref. |T^] found about 



15% violations for an unimproved Wilson action, while results from their 
improved (perfect) action scaled within errors. 

We define the string tension o in terms of a torelon of finite length L by 

a = fji/L + tt/(3L 2 ) , (10) 

and demonstrate scaling in various dimensionless ratios involving a in Ta- 
ble H The extra Coulombic term in Eq. ( |T0"D accounts for long wavelength 
fluctuations, and appears universally in quantized string models when L is 



finite [^3|, Without it, we found that quantities in the last two columns 
for the two different physical lengths independently scale, but show a signif- 
icant discrepancy when compared, providing confirmation for these models. 
Finally, using Mps to determine the scale, we give values for the string ten- 
sion in MeV. Central values from these last two columns are T c f \fa = .649(5) 
and y/a = 467(4) MeV. 



6 Conclusions 

We computed the critical temperature T c for the pure QCD deconfining tran- 
sition on coarse lattices, with spacing a from .33 fm to .15 fm. By comparison 
with charmonium splittings and torelon masses measured at these same spac- 
ings, we demonstrated that a gluon action improved to remove 0(a 2 ) errors 
could measurements nearly free of lattice spacing errors. In particular, that 
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lattices with N s as small as three could give results for T c accurate to 5%. 
This precision is consistent with those found previously for charmonium and 
for the static quark potential. We found central values for T c and the string 
tension a of T c = 300(6) MeV, T c /y/<T = .649(5) and yfd = 467(4) MeV. 
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